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ABSTRACT 

We present the first detailed imaging polarization observations of six SMC and six LMC clusters, 
known to have large populations of B-type stars which exhibit excess Ha emission, to constrain the 
evolutionary status of these stars and hence better establish links between the onset of disk formation 
in classical Be stars and cluster age and/or metallicity. We parameterize the interstellar polarization 
(ISP) along the lines of sight to these twelve clusters, thereby providing a diagnostic of the fundamental 
properties of the dust which characterizes their localized interstellar medium. We determine that the 
ISP associated with the SMC cluster NGC 330 is characterized by a modified Serkowski law with \ max 
^4500A, indicating the presence of smaller than average dust grains. Furthermore, the morphology 
of the ISP associated with the LMC cluster NGC 2100 suggests that its interstellar environment is 
characterized by a complex magnetic field. 

Removing this interstellar polarization component from our data isolates the presence of any intrinsic 
polarization; the wavelength dependence of this intrinsic polarization provides a diagnostic of the dom- 
inant and any secondary polarigenic agents present, enabling us to discriminate pure gas disk systems, 
i.e. classical Be stars, from composite gas plus dust disk systems, i.e. Herbig Ae/Be or B[e] stars. Our 
intrinsic polarization results, along with available near-IR color information, strongly supports the sug- 
gestion of Wisniewski et al. that classical Be stars are present in clusters of age 5-8 Myr, and contradict 
assertions that the Be phenomenon only develops in the second half of a B star's main sequence lifetime, 
i.e. no earlier than 10 Myr. Our data imply that a significant number of B-type stars must emerge onto 
the zero-age-main-sequence rotating at near-critical rotation rates, although we can not rule out the pos- 
sibility that these data instead reveal the presence of a sub-group of the Be phenomenon characterized 
by sub-critically rotating objects. 

Comparing the polarimetric properties of our dataset to a similar survey of Galactic classical Be 
stars, we find that the prevalence of polarimetric Balmcr jump signatures decreases with metallicity. 
We speculate that these results might indicate that either it is more difficult to form large disk systems 
in low metallicity environments, or that the average disk temperature is higher in these low metallicity 
environments. We have characterized the polarimetric signatures of all candidate Be stars in our data 
sample and find ^25% are unlikely to arise from true classical Be star-disk systems. This detection of 
such a substantial number "contaminants" suggests one should proceed with caution when attempting 
to determine the role of evolutionary age and/or metallicity in the Be phenomenon purely via 2-CD 
results. 

Subject headings: Magellanic Clouds — stars: emission-line, Be — circumstellar matter — techniques: 
polarimetric — clustersundividual (Bruck 60, NGC 330, NGC 346, NGC 371, NGC 
456, NGC 458, LH 72, NGC 1818, NGC 1858, NGC 1948, NGC 2004, NGC 2100) - 
starsdndividual (pi-Aquarii, BD+61dl54, MWC349A) 



1. INTRODUCTION 

While the rapid r otation (v eo / V crit ~ 70-80% of 
their critical velocity; iPorteri Il996t iPorter fc Riviniusl 
2003) of classical Be stars has long been specu- 
lated to be the fundamental source driving the pro- 
duction of their geometrically thin circumstellar disks 
dStruvd [l93lt IPorter fc Riviniusl l2003h. recent photomet 



I Grebel fc Chull2000t IKeller et al.|[2000t lOlsen et~all [200ll 
iMcSwain fc Giesl l2005t IWisniewski fc Biorkmanl l2006h 
have suggested secondary mechanisms might contribute 
to the observe d phen o menon. Specifically, iMermi lliod 
(11982ft ; IGrebell (| 19971 ): iFabregat fc Torreionl (|2000D and 
iKellerl ( 20041 ) found that the frequency of the Be 



ric surveys lFeasdl!972t IGrebel. Richtler. fc de B~oerl | 



IGrebell Il997t iDieball fc Grebell 119981 : IKeller et all 



1992; 



1999: 



phenomenon seems to peak in clusters with a main 
sequence turn-off of B1-B2, leading to the sugges- 
tion that the Be phenomenon is enhanced with evo- 
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lutionary age. Several recent o bservational studies 
(iGrebel Richtler. fc de Boerl Il992t iMazzali et all Il996t 
Grebell 119971: iMaeder. Grebel. fc Mermilliod 119991: iKellerl 



2004) have also suggested that the Be phenomenon may 



the true nature of candidate Be stars in six LMC and six 
SMC clusters which were initially identified via photomet- 
ric 2-CDs. In Section 2, we outline our observations and 
data reduction techniques. In Section 3, we detail how we 
identified and removed the interstellar polarization compo- 
nent associated with each of our lines of sight. We develop 
a classification system to characterize the intrinsic polar- 
ization observed in our candidate Be star sample in Section 
4 and also discuss the results of individual clusters. We 
discuss the implications of these results in Section 5 and 
summarize our work in Section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

The imaging polarimetry data presented in this study 
were obtained at the Cerro Tololo Inter-American Ob- 
servatory (CTIOjQ 1.5m. The f/7.5 secondary configu- 
ration was used during our 2001 observing run, which 
yielded a 15 x 15 arc-minute field of view and a 0.44 
arc-second pixel -1 scale, while the f/13.5 secondary con- 
figuration, which provided a 8 x 8 arc-minute field of 
view with a 0.24 arc-second pixel -1 scale, was used dur- 
ing our 2002 observing run. Data were recorded with 
CTIO's standard Cassegrain focus CCD (CFCCD), a 2048 
x 2048 multi-amplifier CCD. The standard CFCCD con- 
figuration was modified by the addition of a rotatable 
half-wave plate followed by a fixed analyzer, the latter of 
which was placed in the upper of the 1.5m's two filter 
wheels. The polymer achromatic wave-plate was manu- 
factured by Meadowlark Optics and provided polarization 
modulation for each 90° rotation of the wave-plate. Un- 
less otherwise noted, each of our fields of view was ob- 
served with the wave-plate rotated through 8 positions 
22.5° apart. The fixed analyzer was a double calcite block 
whose optical axes had been crossed to minimize astig- 
matism and color effects . Further details abou t this po- 
larime t er can be found in |Magalhae s et al.l dl996f). Perevra 
(|2000|) . iMelgareio et~aTl (|200lD . fc iPerevra fc Magalhaes 
(2002J). We placed CTIO's Tek U, B, V, R, and I filters in 
the second filter wheel. A summary of the science clusters 
we observed is presented in Table [T] 

Nightly observations of polarization standard stars were 
obtained to calibrate the polarimetric efficiency of each 
filter and determine the zero-point of the polarization po- 
sition angle in each filter, while nightly observations of un- 
polarized standard stars were obtained to characterize any 
instrumental polarization present. Our unpolarized stan- 
dard data indicate that the polarimeter had an instrumen- 
tal polarization consistent with zero for our 2001 data. The 
B, V, R, and I filters of our 2002 data a lso had an instru- 
ment al polarization consistent with zero (jWisniewski et all 
2003|), while the U filter data had an instrumental polar- 
ization of ~0.1%, which was removed from the data. 

The reduction of these data began with standard image 
processing, including bias and flat field corrections, us- 
ing standard IRAEjf] techniques. Following this initial pro- 
cessing, aperture photometry was performed for all point 
sources in each of the 8 images of an observation set, us- 
ing 10 apertures of size 3-12 pixels for the f/7.5 data and 

7 The Cerro Tololo Inter- American Observatory is operated by the Association of Universities for Research in Astronomy, under contract with 
the National Science Foundation. 

8 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in 
Astronomy, Inc., under contract with the National Science Foundation. 



be more prevalent in low metallicity environments, based 
on comparisons of the fractional Be populations of Galac- 
tic, L MC, and SMC clusters. 

In IWisniewski fc Biorkmanl (|2006l ) (hereafter WB06), 
we used a simple 2-color diagram (2-CD) technique to 
identify the fractional candidate Be population of nu- 
merous Large Magellanic Cloud (LMC), Small Magellanic 
Cloud (SMC), and Galactic clusters in an effort to improve 
the statistical database which has been used to link classi- 
cal Be disk formation with evolutionary age and/or metal- 
licity. WB06 found evidence that the Be phenomenon 
deve lops much earlier than pre viously predicted by the- 
ory (Fa bregat fc Torreionl 120001 ) . i.e before the mid-point 
main sequence lifetime; furthermore, they found evidence 
of an additional enhancement in the fractional Be con- 
tent of clusters with evolutionary age. The increased 
statistics offered by this work, while confirming the previ- 
ously suggested trend of an enhancement in the Be phe- 
nomenon in low metallicity environments, lowered the av- 
erage fractional Be content of SMC clusters from 39% 
(|Maeder. Grebel. fc Mermilliodlll999l) to 32%. 

While the 2-CD has been widely used to link evolution- 
ary age and/or metallicity with Be disk formation, it is 
inherently unclear whether all B-type objects identified as 
excess Ha emitters, i.e. "Be stars" , are truly classical Be 
stars. It has been noted that other B-type objects, such 
as Herbig Ae/Be stars, post-main-sequence B[e] stars, and 
supergiants, may also exhibit Ha emission hence "contam- 
inate" these claimed detections (WB06). 

Polarimetry is a tool which has long been 
used to investigate the circumstcllar environ- 



ments of Be stars (Coyne & Kruszcwski 


1969: 


Poeckert & Marlborough! 11976: 


Covne & McLean 


1982: 


Biorkmanl 11994 McDavidl 1994: 


Quirrcnbach et al. 


1997: 


Wood. Biorkman. & Biorkman '. 


L997; Clarke & Biorkman 


199S|). Electron scattering in the innermost region of clas- 



sical Be circumstellar disks will polarize a small fraction 
of stellar photons; furthermore, pre- and post-scattering 
attenuation by hydrogen atoms within the inner disk re- 
gion will imprint a characteristic wavelength dependence 
onto this polari zation signature if the inn er disk is suf- 
ficiently dense (|Wood fc Biorkmanl |1995[ ). The wave- 
length dependence of intrinsic polarization originating 
from classical Be stars, charact erized by geometrically 
thin, gaseous circumstellar disks. (lOuirrenbach et al.ll997t 
Wood, B iorkman. fc Biorkmanl [19971 : IClarke fc B iorkman 



1998) is significantly different than that expected from 



the dustier circumstellar environments of Herbig Ae/Be 
stars, post-main - sequence B\e] sta rs, and supergiants 
(|Vrba et al.lll979l : iMever et al.ll2002t ). Thus for non-pole- 
on disk geometries, polarimetry can be employed to in- 
vestigate circumstellar environments and discriminate the 
evolutionary nature of sources. 

In this paper, we use imaging polarimetry to investigate 
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5-14 pixels for the f/13.5 data. The linear polarization of 
each source was extracted from a least squares solution 
of the difference amplitudes in the 8 wave- plate positions 
(A), using the PCCDPACK software suite (|Perevrall2000l: 
iPerevra fe Magalhaesll2002ft . Both the expected photon 
noise errors and the actual measured errors, defined as the 
residuals at each wave-plate position with respect to the 
expected cosiA curve, were also calculated. With few 
exceptions, these two errors were consistent with one an- 
other. 

3. INTERSTELLAR POLARIZATION 

The polarization observed in our raw data is comprised 
of the superposition of an interstellar component, at- 
tributable to the dichroic absorption of starlight by par- 
tially aligned dust grains along each line of sight, and pos- 
sibly an intrinsic component, attributable to the asymmet- 
ric illumination of circumstcllar material and/or the illu- 
mination of an asymmetrical distribution of circumstellar 
material. Characterizing and isolating each of these com- 
ponents, although often technically challenging, yields two 
unique diagnostics. Isolating the interstellar polarization 
component can facilitate detailed studies of the local in- 
terstellar medium (ISM) and provide constraints on fun- 
damental properties such as grain size distributions and 
shapes, as well as diagnosing the magnitude and direction 
of local magnetic fields. Isolating the intrinsic polarization 
component enables investigations of the distribution and 
chemistry of scatterers which comprise unresolved circum- 
stellar environments. We stress that characterizing the 
polarizing agent (i.e. gas or dust) responsible for produc- 
ing an intrinsic polarization signal is not possible simply 
via inspection of the wavelength dependence of the vector 
sum of interstellar plus intrinsic (i.e. the total) polariza- 
tion data. 

Fortunately, our efforts to identify and separate the in- 
terstellar and intrinsic polarization components of our data 
are simplified by the fact that we have observed rich stellar 
associations. Most of the stars in these clusters should be 
normal main sequence stars lacking any type of asymmet- 
rical circumstellar envelope, hence they should not exhibit 
an intrinsic polarization component. Furthermore, we can 
assume that all of our stars are located at the same gen- 
eral distance, since they are members of a cluster popula- 
tion, and that the properties of the interstellar medium do 
not significantly vary over the angular extent of our clus- 
ters. Hence each cluster in our sample provides us with a 
multitude of suitable "field stars" whose polarization we 
can simply average to estimate the interstellar polarization 
(ISP) along eac h line of sight, via the co mmonly used field 
star technique (jMcLean fc Brownlll978l ). 

In Figure [T] we show the total V-band polarization of all 
sources in the LMC cluster NGC 1818 having polarimetric 
signal-to-noise ratios (p/c P ) greater than 5.0. The distinc- 
tive grouping of most objects in the Stokes Q-U diagram 
of this figure (panel a) and the narrow polarization (panel 
d) and position angle (panel c) histograms demonstrates 
that most of these objects do lack significant intrinsic po- 
larization components and only have an ISP component. 
After excluding objects outside of the dominant trends 
found in these histograms with an iterative method, we 
calculated the weighted average and standard deviation of 



these data to determine a preliminary estimate of this clus- 
ter's V-band ISP, as shown in Figure [2] The wavelength 
dependence of t hese estimates clearly follow the empirical 
Serkowski law ([Serkowski et al.lll975ft commonly used to 
parameterize an ISP . Over-plotted in Figure [H is a mod- 
ified Serkowski law (| Serkowski et al.l 119751 : IWilking et alj 
11982ft which we deemed best represented the observational 
data. The polarization at the midpoint wavelength loca- 
tion of the U, B, V, R, and I filters was then extracted from 
this curve and served as the final estimate of the ISP. The 
modified Serkowski parameters determined for all of our 
lines of sight, along with the extracted filter ISP values, 
are summarized in Table [5J 

The ISP values compiled in Table [2] are relevant to our 
efforts to isolate the intrinsic polarization components as- 
sociated with our target stars; however, they do not yet 
yield any direct information regarding the properties of 
the interstellar dust grains which reside in their parent 
LMC and SMC clusters, as these ISP values are com- 
prised of both foreground Galactic ISP and Magellanic 
Cl oud ISP contri b ution s. Fo llowing the techniques used 
by IClavton etai] (jl983l ) and iRodrigues et alj (119971), w e 
used the interstellar polarization maps of lSchmidtl l| 1976ft . 
to identify and remove this foreground Galactic ISP con- 
tribution (see Table [3]), hence determine estimates of the 
ISP intrinsic to our LMC and SMC clu sters (see Table 
U). We used a modified Serkowski-law (| Serkowski et alj 
1975; Wilkin g~et al.lll982ft to extrapolate these maps over 



the wavelength range of our dataset, and assumed a nom- 
inal Galactic \ m ax value of 5500A to characterize the pa- 
ramet er K. The erro rs cited in the Galactic polarization 
maps (|SchmidtlH976ft and those present in our total ISP 
estimates (Table d]) were propagated to produce the error 
estimates cited in Table SJ 

Most of the field stars used to derive our total ISP es- 
timates were not highly reddened objects, hence the in- 
trinsic LMC and SMC ISP values listed in Table H all are 
characterized by modest polarization amplitudes, which 
are often on the order of the errors of our data. In spite 
of this fact, trends in these data are clearly present. As 
expected, the magnitude of UBVRI polarization towards 
each cluster follows a Serkowski-like wavelength depen- 
dence at a wavelength independent position angle. This 
lack of position angle rotation indicates that our Galac- 
tic interstellar polarization correction was reasonable and 
that we are looking at a single magnetic field orientation 
towards each of our clusters. All LMC clusters exhibit 
similar ISP properties, with polarization magnitudes rang- 
ing from ^0.2-0.5% at a position angle of ^25-45 degrees; 
similarly, all SMC clusters also exhibit clear evidence of 
sharing common ISP properties, with polarization magni- 
tudes ranging from ^0.3-0.6% at a position angle of ^120- 
150 degrees. The shallow curvature of the Serkowski-law 
dependence of these data, in combination with the mod- 
erate level of uncertainty present, make identifying sys- 
tematic differences in the interstellar dust grain properties 
of these clusters difficult. The total ISP (Table ^ and 
SMC ISP (Table gj) of the cluster NGC 330 does how- 
ever show suggestive evidence of being characterized by 
a short wavelength X max value of ^4500 A (Figure [3]). 
Such an effect is commonly attributed to the presence of 
small dust grains, and has been previously observed in 
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other SMC sight-lines (jRodrigues et al.lll997t ). In a fu- 
ture paper, we will examine the wavelength dependence 
of moderately and highly reddened objects in our dataset 
to perform a more detailed investigation of the interstellar 
medium properties of these LMC and SMC clusters, and 
compliment earlier studi es which parameterize d the ISM 
properties of other S MC (jRodrigues et al.lll997l) and LMC 
(| Clayton et~ailll983D sight-lines. 

3.1. NGC 2100 

The total polarization vectors of our observations of the 
LMC cluster NGC 2100 clearly exhibit evidence of a col- 
lective, complex morphology (see Figure 2]). While Figure 
H] only presents I-band data, all of our other filters exhibit 
similar alignment patterns. Recall that we expect most 
cluster objects should not exhibit an intrinsic polarization 
component; thus we suggest that these systematic mor- 
phological changes are related to changes in the magnetic 
field properties within our field of view. A substantial dis- 
cussion of t he magnetic field p ropert ies of NGC 2100 is 
presented in IWisniewski et all (|2007aD . For the purposes 
of the present study, we will only discuss our efforts to pa- 
rameterize and remove the ISP along the line of sight to 
this cluster. 

We identified three spatial regions in our field of view, 
corresponding to the region around the cluster core, the 
region to the south of the cluster core, and the region to 
the north of the core, which displayed unique ISP char- 
acteristics. In Figures [5J H3 and [7| we show the B-band 
polarization vector maps of what we have defined as re- 
gions 1, 2, and 3 respectively in the total ISP along the 
line of sight to NGC 2100, overlaid on DSS-2 blue images. 
Having identified these distinct regions, we then extracted 
total ISP estimates following the technique previously de- 
scribed. We calculated a er~ 2 weighted average of all ob- 
jects in each filter served as an initial ISP estimate, deter- 
mined the modified Serkowski-law parameters which best 
represented the data, as seen in Figure [5] for region 1, and 
extracted final ISP values (Table [5} for each region from 
these Serkowski curves. The locations of the candidate Be 
stars in NGC 2100 were then correlated to these 3 ISP re- 
gions to determine the ISP correction each should receive. 
We found NGC 2100:KWBBe 102, 630, and 712 resided in 
area 2, NGC 2100:KWBBe 535, 797, and 1033 resided in 
area 3, and the rest of the candidate Be stars resided in 
area 1. 

4. INTRINSIC POLARIZATION 

We subtracted the ISP from the observed polarization 
to isolate any intrinsic polarization components present. 
The location of candidate Be stars in our fields of view 
were identified by careful correlation with literature co- 
ordinates and finder charts. Using several PCCDPACK 
routines, we meticulously examined the extracted polar- 
ization for each of the candidate Be stars in our sample 
to search for undcsircd contamination by a) the ordinary 
or extraordinary images of nearby objects and b) cosmic 
ray hits or uncorrected bad pixels. The total polarization 
of all candidate Be stars which did not suffer from these 
contamination issues are tabulated in Table O while the 
intrinsic polarization of these objects are tabulated in Ta- 
ble[7] 



As discussed in the Introduction, electron scattering in 
the pure gas disks of classical Be stars will polarize a small 
fraction of stellar photons, producing a well known wave- 
length dependent intrinsic polarization signature. We em- 
phasize that as polarization is a vector quantity, identi- 
fying these signatures a) is only possible via inspection 
of intrinsic polarization data; and b) requires either the 
simultaneous inspection of both the wavelength depen- 
dence of the polarization magnitude and position angle, 
or alternatively inspection of the behavior of the intrin- 
sic data on a Stokes Q-U diagram. An observed wave- 
length independent intrinsic polarization magnitude and 
position angle across the entire UBVRI wavelength regime 
is the expected signature of pure electron scattering, such 
as that expected from a gaseous classical Be disk. More- 
over, pre- or post-scattering absorption will superimpose 
the wavelength dependent signature of hydrogen opacity 
on th is signal, creatin g a "saw-tooth" polarization signa- 
ture (|Biorkmanll2000l ). if enough absorption events occur. 
Specifically, this "saw-tooth" signature includes the pres- 
ence of abrupt jumps in the magnitude of intrinsic po- 
larization at the Balmer and Paschen limits, as well as a 
wavelength independent intrinsic polarization position an- 
gle across the entire UBVRI wavelength regime. Examples 
of these two signatures of classical Be circumstellar disks 
are given in Figure [9l which presents multi-epoch observa- 
tions of the known classical Be star pi Aquarii as observed 
by the University of Wisconsin's HPOL spectropolarime- 
ter. 

We have contemporaneously analyzed the Stokes Q-U 
diagrams of all candidate Be stars' total and intrinsic po- 
larization components to search for evidence of these sig- 
natures. We have developed a conservative 4-point classifi- 
cation scale to rank the likelihood that candidate Be stars 
are truly classical Be stars. Results of this classification 
are summarized in Table [5] for individual objects and in 
Table [9] for the net results of entire clusters. The design of 
our classification system emphasized the identification of 
objects exhibiting a "saw-tooth" intrinsic polarization sig- 
nature and a pure electron scattering signature across the 
entire UBVRI wavelength regime. Specifically, we defined 
our classification system by: 

• Type-1: Objects which are most likely classical Be 
stars; 

• Type-2: Objects whose polarimetric properties are 
not inconsistent with those expected from classical 
Be stars; 

• Type-3: Objects which are unlikely to be classical 
Be stars; and 

• Type-4: Objects which are highly unlikely to be 
classical Be stars. 

Stars in our sample which displayed, to within 3 c, 
a "sawtooth-like" polarization Balmer jump (BJ) along 
with a wavelength independent polarization position an- 
gle across the entire range of available data (see Figure [TU]) 
were assigned a designation of type-1. We claim that all 
type-1 objects are most likely classical Be stars. Note that 
we offer a detailed discussion of other objects which could 
exhibit somewhat similar behavior over part of the UB- 
VRI wavelength regime in Section 15.3.11 However, as we 
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discuss in Section [5.3.31 analysis of ancillary data demon- 
strates that it is highly dubious that stars classified as 
Type-1 polarimetric sources are anything other than clas- 
sical Be stars. 

We assigned a designation of type-2 to objects whose 
polarimetric properties were not inconsistent with that ex- 
pected from a classical Be star-disk system. As previously 
discussed, classical Be disks of sufficiently low density will 
not leave an imprint of hydrogen opacity in their intrin- 
sic polarization signals; hence, stars which exhibited, to 
within 3cr, a wavelength independent electron scattering 
(ES) polarization signature (see Figure ITT]) received a des- 
ignation of type-2. A small number of objects exhibited 
a nearly wavelength independent polarization magnitude 
along with a minor wavelength dependence in their po- 
larization position angles; we designated these objects as 
type-2 as we believe such stars most likely are exhibit- 
ing ES signatures modified by a slight under- or over- 
correction of their ISP components. We can not rule 
out that this minor wavelength dependence might also 
be produced by the additional presence of an optically 
thin dust disk, as present in post-main-s equence B[e] stars 
(|Magalhaesl [l992 ; Melgareio et alj|2001h . However, as we 
show in Section 15.3.31 analysis of ancillary data demon- 
strate that most type-2 objects which we claim exhibit 
electron-scattering (ES) signatures are likely to be clas- 
sical Be stars, and substantially less likely to be systems 
characterized by composite gas plus dust disks (i.e. Herbig 
Ae/Be, B[e] stars). Finally, we assigned all objects which 
appeared to be intrinsically unpolarized (< 0.3% polar- 
ization), to within 3er, a type-2 designation as this is the 
expected signature of pole-on or nearly pole-on classical Be 
stars. We caution the reader that we are not 100% certain 
that these unpolarized objects are classical Be stars as: a) 
stars without gaseous disks, spuriously detected as excess 
Ha emitters on 2-CDs, will also exhibit zero net intrinsic 
polarization and b) the noise present in the observations 
of fainter targets may preclude us from clearly identifying 
them as "contaminants" . 

We assigned the designation of type-3 to type-4 to ob- 
jects whose polarization, to within 3c, appeared to be 
inconsistent with the aforementioned signatures expected 
from a classical Be star-disk system (see e.g. Figure [T2l. 
We briefly discuss some of the major types of atypical po- 
larimetric signatures we observed. 

In Table [SJ we assigned a designation of type 3-4 to sev- 
eral objects whose intrinsic polarization exhibited signs of 
a 90° position angle reversal (Figure [TT]). which is a signa- 
ture of a dusty bipolar nebu l a geometry (see Sect ion [5JLU 
see also lSchmidt et al.lll992tlSchulte-Ladbeck et al.l fl992). 
From an inspection of the Stokes Q-U diagram of such an 
object (Figure [TT]) . it is clear that even if our initial ISP 
correction was grossly miscalculated, the wavelength de- 
pendent polarization would still be inconsistent with that 
expected from classical Be stars. While it would be inter- 
esting to further probe the circumstellar environments of 
such objects with follow-up investigations, for the purposes 
of this paper we merely remark that they are unlikely to 
be classical Be stars. 

Several of our candidate Be stars exhibited very large (1- 
3.5%), complex intrinsic polarization signatures (see Fig- 
urc !13[) . Given the steep drop in their polarization magni- 



tude at short optical wavelengths, along with subtle indica- 
tions of a corresponding position angle rotation, it is pos- 
sible that these objects are characterized by dusty bipolar 
nebulae whose polarization signature exhibits a position 
angle reversal at UV wavele ngths, similar to HD 45677 
( Schulte-Lad beck et al.|[l~992] ). Alternatively, the polariza- 
tion of these objects could be interpreted as following a 
Serkowski-like wavelength dependence characterized by a 
large X m ax value. As our data have already been cor- 
rected for the average ISP associated with each of our 
clusters, this latter interpretation would require these ob- 
jects to be situated in region of patchy dust, likely popu- 
lated by larger grains given the lon g wavelength values of 
A mnx , i.e. > 7 000 A in Figure [T3l (IRodrigues et alJll997t 
Whi ttet et al.l[l99l . iSerkowski et al.l(|1975h described the 
expected magnitude of interstellar polarization by the for- 
mulation 3Eb~v < Pmax < QEb-v] hence, an additional 
reddening of 0.3 - 1.0 E(B — V would be required to pro- 
duce the measured p ma x of these objects, ~3%. The 
observed colors of at least some of these objects (NGC 
1948:KWBBe 98, V-I = 0.67; Keller et al. 1999) is likely 
sufficient to produce the amount of purported additional 
ISP, although other objects (NGC 1948:KWBBe 246, V-I 
= 0.05; Keller et al. 1999) clearly do not meet this criteria. 
While follow-up optical or infrared spectroscopy of these 
anomalous candidates would help to determine whether 
they exhibit a dusty bipolar nebula geometry or are lo- 
cated in a region of patchy dust, for the purpose of this 
paper we merely stress that it is clear that these objects 
are not likely to be classical Be stars. 

We now offer discussion of the polarimetric properties 
of candidate Be stars in individual clusters. Recall that 
all candidates which exhibited evidence of contamination, 
predominantly from nearby neighbors, have been excluded 
from our analysis. We further note that the effective lim- 
iting magnitude of our polarimetric dataset was lower (i.e. 
brighter) than that considered in WB06. 

4.1. LMC Clusters 

4.1.1. NGC 1818 

iKeller et al] (|l999t ) identified 40 candidate Be stars asso- 
ciated with the LMC cluster NGC 1818 and its surround- 
ing field. The cluster is densely populated, and because of 
image overlap issues, we were only able to retrieve polari- 
metric information for 18 of these candidate Be stars. 4 of 
these 18 (22%) showed intrinsic polarization Balmer jump 
(BJ) signatures, NGC 1818:KWBBe 69, 82, 137, 243. We 
classified 12 of the 18 candidates (67%) as type-2 objects, 
and found 5 of the type-2 objects (i.e. 5 of the 18 stars, 
28%, in the total population) exhibited clear evidence of 
an electron scattering (ES) polarization signature. 2 of the 
18 candidates (11%), NGC 1818:KWBBe 47 and 381, were 
deemed unlikely to be classical Be stars. While the B, V, 
R, and I filter polarization of NGC 1818:KWBBe 47 were 
consistent with an ES signature, the U filter exhibited a 
significant position angle rotation which was inconsistent 
with an ES origin. Although noisy, the observation of 
NGC 1818:KWBBe 381 also did not follow an ES wave- 
length dependence, prompting us to assign it an "unlikely 
Be star" designation. 

4.1.2. NGC 1948 
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iKeller et al.l (|l999f ) identified 27 candidate Be stars as- 
sociated with the cluster NGC 1948 and its nearby field. 
We were able to extract polarimetric information on 22 of 
these candidates from our data set. 6 of these 22 candi- 
dates (27%) exhibited a polarization Balmer jump, NGC 
1948:KWBBe 71, 75, 102, 153, 172, and 240. We classi- 
fied 12 of the 22 candidates (55%) as type-2 objects, and 
found 6 of the type-2 objects (i.e. 6 of the 22, 27%, in 
the total population) exhibited a definite electron scat- 
tering polarization signature, NGC 1948:KWBBe 62, 92, 
101, 157, 326, and 790. We note that NGC 1948:KWBBe 
92 does show a hint of a small polarization Balmer jump; 
however, we opted to describe its polarization more con- 
servatively, i.e. as having an electron scattering signature. 
We found 4 of the 22 stars (18%) had polarimetric prop- 
erties which suggested that they were unlikely to be clas- 
sical Be stars. As previously discussed, 3 of these stars, 
NGC 1948:KWBBe 98, 183, and 246 exhibited intrinsic 
polarization signals which seemed to follow a Serkowski- 
like wavelength dependence. NGC 1948:KWBBe 91 was 
also deemed unlikely to be a classical Be star due to the 
wavelength dependent nature of its polarization position 
angle. 

4.1.3. NGC 2004 

IKeller et~ al. (1999) identified 67 candidate Be stars asso- 
ciated with the LMC cluster NGC 2004 and its surround- 
ing field. We extracted polarimetric information for 43 of 
these candidates, and found that 9 of the 43 (21%), namely 
NGC 2004:KWBBe50, 87, 96, 103, 106, 152, 211, 347, 377, 
showed an intrinsic polarization Balmer jump. 28 of the 
43 candidates (65%) were classified as type-2 objects, and 
note that 10 of these 28 type-2 objects (i.e. 10 of the 43, 
23%, in the total population), NGC 2004:KWBBe 91, 203, 
276, 323, 344, 441, 624, 717, 1175, and 1421 exhibited an 
electron scattering polarization signature. We found 6 of 
the 43 (14%) candidates in NGC 2004 were unlikely to 
be classical Be stars based on their intrinsic polarization. 
Of these objects, we note that the 3 filters of polarization 
extracted for NGC 2004:KWBBe 1315 show suggestive ev- 
idence of a polarization position angle flip, a feature which 
is not expected in classical Be star-disk systems. 

4.1.4. LH 72 

WB06 identified 50 candidate Be stars in the LMC clus- 
ter LH 72, and designated 11 of these detections as tenta- 
tive. We were able to obtain polarimetric information for 
34 of these stars: 1 of the 34 (3%), LH 72WBBe 5, showed 
a polarization Balmer jump and is most likely a bona-fide 
classical Be star. We classified 22 of the 34 (65%) stars as 
type-2 objects, and found 5 of these 22 type-2 objects (i.e. 
5 of the 34, 15%, of the total population), LH 72WBBc 
13, 15, 26, 27, and 33 exhibited an electron scattering po- 
larization signature. We suggest that 11 of the 34 (32%) 
of candidates in LH 72 are unlikely to be classical Be stars 
based on their intrinsic polarization properties. Of these 
unlikely Be stars, we found that the intrinsic polarization 
of LH 72WBBc 9 followed a Serkowski-like wavelength 
dependence, with a X ma x value > 6000 A. We extracted 
polarimetric information for 7 of the 11 stars designated 
as possible candidate Be stars by WB06. 5 of these 7 stars 
were classified as type-2 objects, while 2 of the 7 appear 



unlikely to be classical Be stars based on their intrinsic 
polarization signatures. 

4.1.5. NGC 1858 

WB06 identified 39 candidate Be stars in the LMC clus- 
ter NGC 1858, and we were able to extract polarimetric in- 
formation for 27 of these 39 stars. 3 of the 27 (11%), NGC 
1858WBBe 6, 9, and 20, exhibited polarization Balmer 
jumps. We classified 13 of the 27 (48%) as type-2 objects, 
and found 2 of these 13 type-2 stars (i.e. 2 of the 27, 7%, 
of the total population) exhibited an electron scattering 
polarization signature. 11 of the 27 (41%) candidates in 
NGC 1858 appear unlikely to be classical Be stars based on 
their polarimetric signatures. As previously discussed, 2 of 
these unlikely Be stars, NGC 1858 WBBe 3 and 12 had an 
intrinsic polarization which follow a Serkowski-like wave- 
length dependence, characterized by a long X max value. 4 
of the 39 photometrically identified candidate Be stars in 
NGC 1858 were judged to be possible detections in WB06: 
we were able to detect 3 of these 4 polarimetrically. We 
classified 2 of the 3 to be type-2 objects, and the remaining 
object, NGC 1858:WBBe 9, as a type-1 object. 

4.1.6. NGC 2100 

IKeller et al.l (|1999t) identified 61 candidate Be stars as- 
sociated with the LMC cluster NGC 2100. We were able 
to extract polarimetric information for 35 of the 61 candi- 
dates (57%). 8 of the 35 stars (23%), NGC 2100:KWBBe 
79, 97, 436, 619, 635, 705, 770, and 797, exhibited po- 
larization Balmer jumps and hence are most likely to be 
classical Be stars. 20 of the 35 stars (57%) were classified 
as type-2 stars, and we found 16 of these 20 type-2 stars 
(i.e. 16 of the 35, 46%, of the total population) exhibited 
an electron scattering intrinsic polarization signature. We 
suggest that 7 of the 35 (20%) candidates in NGC 2100 are 
unlikely to be classical Be stars based on their observed 
intrinsic polarization signatures. Note that the intrinsic 
polarization of NGC 2100:KWBBe 321 seems to follow a 
Serkowski-like wavelength dependence, while the unique 
intrinsic polarization signatures of NGC 2100:KWBBe 111 
and 219 (Figure fT4|) suggests these objects might be dust- 
disk systems. 

4.2. SMC Clusters 

4.2.1. NGC 346 

IKeller et al.l (|1999T ) identified 48 candidate Be stars in 
the vicinity of the SMC cluster NGC 346: we have ob- 
tained polarimetric information for 33 of these objects. 
8 of the 33 objects (24%), NGC 346:WBBe 85, 93, 191, 
236, 374, 445, 468, and 529, exhibited polarization Balmer 
jumps and hence are most likely classical Be stars. We 
classified 22 of the 33 stars (67%) as type-2 stars, and 
found 7 of these 22 type-2 objects (i.e. 7 of the 33, 21%, 
of the total population) exhibited a clear electron scat- 
tering intrinsic polarization signature. We suggest that 
3 of the 33 candidates in NGC 346 (9%) are unlikely to 
be classical Be stars based on their polarimetric signa- 
tures. Although not detected via our polarimetric survey, 
we note that NGC 346:KWBBe 13 and NGC 346:KWBBe 
200 should not be considered to be classical Be stars; the 
former (NGC 346:KWBBe 13) is the well-known Wolf- 
Rayet/LBV HD 5980 and the later (NGC 346:KWBBe 
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200) has r ecently been shown to be the fifth known B[e] in 
the SMC (|Wisniewski et al.ll2007bh . 

4.2.2. NGC 371 

WB06 identified 129 candidate Be stars in NGC 371. 
We have obtained polarimetric information for 73 of these 
targets, and found 10 of these 73 stars (14%), NGC 
371:WBBe 2, 3, 4, 5, 6, 10, 13, 18, 24, and 31 exhibited 
polarization Balmcr jumps, indicating they are most likely 
classical Be stars. We classified 49 of the 73 stars (67%) as 
type-2 objects, and remark that 6 of these 49 type-2 stars 
(i.e. 6 of the 73, 8%, of the total population) exhibited a 
clear electron scattering polarization signature. Note that 
many of the fainter candidate Be stars included in our 74 
detections were not observed at high signal-to-noise lev- 
els; thus for most of these objects we can only say that 
to within 3(7, their polarimetric properties are not incon- 
sistent with those of classical Be stars. We suggest that 
14 of the 73 (19%) candidates in NGC 371 are unlikely to 
be classical Be stars based on their polarimetric proper- 
ties. WB06 suggested that 11 of the 130 photometrically 
identified candidates should be viewed as "possible detec- 
tions" ; however, we were only able to extract polarimetric 
information for 2 of these 11 objects, NGC 371:WBBe 64 
and 87. We classified both as type-2 stars. 

4.2.3. Bruck 60 

We were able to extract polarimetric information for 18 
of the 60 candidate Be stars identified photometrically by 
WB06. We found 1 of the 18 (6%), Bruck 60WBBe 6, 
showed a polarization Balmer jump, while we classified 11 
of the 18 stars (61%) as type-2 objects. 8 of these 11 type- 
2 stars (i.e. 8 of the 18, 44%, of the total population) 
exhibited an electron scattering intrinsic polarization sig- 
nature. We suggest that 6 of the 18 (33%) candidates 
in Bruck 60 are unlikely to be classical Be stars based on 
their polarimetric signatures. WB06 noted that 5 of the 26 
photometrically identified candidate Be stars in Bruck 60 
should be considered "possible detections" . We were only 
able to detect 1 of these 5 stars polarimetrically, Bruck 
60WBBe 21, and classified it as a type-2 object which 
exhibited clear evidence of an electron scattering polariza- 
tion signature. 

4.2.4. NGC 456 

23 candidate Be stars in NGC 456 were identified in the 
photometric survey of WB06, and we were able to extract 
polarimetric information for 14 of these 23 candidates. Al- 
though we did not observe any objects with polarimetric 
Balmer jumps, we did classify 9 of the 14 stars (64%) as 
type-2 objects. 1 of these 9 type-2 stars (i.e. 1 of the 14, 
7%, of the total population), NGC 456WBBe 1, exhibited 
an electron scattering polarization signature. We suggest 
that 5 of the 14 (36%) candidates in NGC 456 are unlikely 
to be classical Be stars based on their observed intrinsic 
polarization components. WB06 suggested 1 of NGC 456's 
23 photometrically identified candidate Be stars should be 
viewed as a "possible detection" : we marginally detected 
this object, NGC 456WBBe 21, and suggest that it is 
unlikely to be a classical Be star. 



4.2.5. NGC 458 

30 candidate Be stars were photometrically identified in 
NGC 458 by WB06, and we were able to extract polari- 
metric information for 10 of these 30 candidates. None 
of the 10 detected objects exhibited polarimetric Balmer 
jumps or clear evidence of having an electron scattering 
signature. We classified 6 of the 10 stars (60%), and sug- 
gest that 4 of the 10 (40%) are unlikely to be classical 
Be stars. We were able to extract polarimetric informa- 
tion for 1 of the 2 initially identified as "possible detec- 
tions" by WB06, NGC 458:WBBel9. We classified this 
star as a type-2 object; furthermore, we found it exhibited 
suggestive evidence of a small polarimetric Balmer jump, 
although the data were very noisy. 

4.2.6. NGC 330 

iKeller et alj (|1999T ) identified 76 candidate Be stars in 
the vicinity of NGC 330 and we were able to extract po- 
larimetric information for 41 of these objects. While we 
did not observe any of the detected candidates to have 
a polarimetric Balmer jump, we did classify 24 of the 41 
(59%) candidates as type-2 stars, and note that 9 of these 
24 type-2 objects (i.e. 9 of the 41, 22%, of the total pop- 
ulation) exhibited electron scattering polarimetric signa- 
tures. We suggest that 17 of the 41 (41%) candidates are 
unlikely to be classical Be stars based on their intrinsic 
polarimetric signatures. 

5. DISCUSSION 
5.1. Analysis of Intrinsic Polarization Statistics 

It is clear from Table [5] that the percentage of SMC and 
LMC candidate Be stars which exhibit polarization Balmcr 
jumps, electron scattering signatures, or were deemed un- 
likely to be classical Be stars varies significantly from one 
cluster to the next. To better assess the global trends 
of our dataset, we extracted simple statistical mean and 
median values from Table O as summarized in Table 1101 
From Table \TU[ we see that at least 25% of our dataset 
is populated by objects which appear unlikely to be clas- 
sical Be stars. The presence of such a large number of 
likely "contaminants" is interesting, as it confirms one of 
the initial hypotheses of this project: it is dangerous to 
assume, a-priori, that all objects identified as excess Ha 
emitters via 2-CD surveys are classical Be stars. We thus 
suggest caution should be exercised when attempting to 
ascertain the role evolutionary age and/or metallicity play 
in the development of the Be phenomenon purely via the 
analysis of 2-CD data. 

Assuming these objects are truly "contaminants" , we 
recalculated the statistics of our intrinsic polarization 
dataset after removing these objects from consideration. 
The results of this exercise are listed in Table [TTJ The 
median prevalence of polarization Balmer jumps in our 
SMC/LMC dataset is typically between 20-25% (Table 
lll[) . Polarization Balmer jumps appear to be significantly 
less prevalent in our SMC clusters (4%) than in our LMC 
clusters (25%); however, this statistic is strongly influ- 
enced by the null detection of such signatures in three 
SMC clusters. 

Sin ce 1989, the HPOL spectropolarimeter (|Wolff et alj 
Il996f ) mounted on the University of Wisconsin's Pine Bluff 
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Observatory (PBO) has been monitoring the optical spec- 
tropolarimetric properties of a sample of 73 known Galac- 
tic Be stars, and it provides a statistically significant po- 
larimetric dataset which we will use to aid the interpre- 
tation of our results. As yet, most of the ~800 HPOL 
observations of classical Be stars have not had their ISP 
components identified and removed; hence, we can not 
readily compare their intrinsic polarization characteristics 
to our data. However, we can determine the prevalence of 
polarization Balmer jumps in the HPOL database. Data 
ob tained prior to 1995 ha ve been compiled in catalog form 
bv lBiorkman et ail (|2000D . and observations characterized 
by a polarization Balmer jump have already been flagged 
in the catalog. We have similarly analyzed all observations 
obtained from 1995-2004 to identify the presence of polari- 
metric Balmer jumps in these data. We then compiled all 
HPOL observations of individual classical Be stars into two 
categories: a) classical Be stars which had never exhib- 
ited a polarimetric Balmer jump in any PBO observation; 
and b) classical Be stars which had exhibited a polarimet- 
ric Balmer jump in at least one observation. The results 
of this classification procedure are given in Table 111! 

We found 31 of the 73 (42%) observed Galactic Be stars 
exhibited a polarimetric Balmer jump for at least part of 
the time-frame of the HPOL survey. Amongst early-type 
Be stars (09-B5) in the HPOL database, 22 of 45 (49%) 
exhibited polarimetric Balmer jumps. Note that further 
restricting the sample to include only 09-B3 type Be stars 
does not appreciably affect the observed prevalence of po- 
larization Balmer jumps, as 17 of 40 (43%) exhibited such 
a signature. In contrast, the prevalence of polarization 
Balmer jumps does seem to be noticeably lower for later- 
type (B6-A0) Be stars, as we found only 8 of 26 (31%) 
exhibited evidence of a polarimetric Balmer jump. 

The frequency of polarization Balmer jumps appears to 
decrease with metallicity: ^45% of Galactic Be stars in the 
HPOL database, ~25% of LMC Be stars from the present 
study, and ~4% of SMC Be stars from the present study 
exhibit polarization Balmer jumps. We suggest several 
possible explanations for these results. 

1. It is possible that the average disk properties of 
classical Be stars are fundamentally different in the 
low metallicity environments of the SMC and LMC 
as compared to our Galaxy. The presence of intrin- 
sic polarization in classical Be stars, as well as its 
wavelength dependence, is a function of disk incli- 
nation angle, disk density, and the effective tem- 
perature of disk material. Thus the observed lower 
frequency of polarization Balmer jumps in our data 
suggest either that it might be harder to form mas- 
sive disk systems in low metallicity environments or 
that the average disk temperature is higher in these 
lower metallicity environments, hence decreasing 
the amount of pre- and post-scattering absorption 
by neutral H I. 

2. Our polarimetric observations provided a one time 
sampling of LMC and SMC candidate Be stars, 
while the HPOL survey obtained a larger number 
of observations, albeit at quasi-random intervals, of 
a population of known Be stars. Owing to the vari- 
able nature of the Be phenomenon, we do not a 



priori know how often any individual Be star will 
exhibit any of the possible polarization signatures; 
hence, the different sampling errors associated with 
these two databases might explain at least some of 
the observed differences in the polarization Balmer 
jump frequencies. 

3. The statistics presented in Table [TT] were based on 
removing a significant number of objects deemed 
unlikely to be classical Be stars from considera- 
tion, i.e. type-3 and type-4 objects. Recall however 
that our definition of type-2 objects only required 
their polarimetric properties to be not inconsistent 
with those expected from classical Be stars. It is 
possible that a population of "contaminants" , ei- 
ther a) unpolarized objects which aren't pole-on or 
near pole-on classical Be stars; or b) non-classical 
Be stars which were not identified as such due to 
low data quality, still reside within our collection 
of type-2 objects, hence underestimate our polar- 
ization Balmer jump fraction frequency. If such a 
scenario is true, it would further suggest that ex- 
treme caution should be exercised when assessing 
the nature of stellar populations purely via the use 
of photometric 2-CDs. 

5.2. Identification of "Contaminants" in 2-CD Detected 
Candidate Be Stars 

WB06 identified a large number of B-type objects hav- 
ing excess Ha emission via photometric 2-CDs and sug- 
gested that these stars be classified as candidate Be stars. 
D. Gies (2006, personal communication) noted that sev- 
eral candidate Be stars in the WB06 dataset had colors 
redder than that expected from classical Be stars, and sug- 
gested that the frequency of observed blue (or red) excess 
Ha emitters in the 2-color diagrams presented in WB06 
increased with the total number of stars of a given color 
range. Furthermore, he postulated that if these red ex- 
cess emitters were a) a manifestation of extremely young 
stars with their birth disks still intact, b) red- type giants 
or supergiants, or c) false detections owing from the effects 
of bright background H II emission; measuring the preva- 
lence of these excess emitters could serve as an estimate of 
the number of "contaminants" which affect the blue pop- 
ulation of excess Ha emitters, i.e. the candidate Be star 
population. 

We examined the 2-CD photometry of 6 clusters, com- 
mon to both this survey and that of WB06 and spanning 
a wide range of ages (very young, young, old; WB06) and 
metallicities (SMC, LMC), and determined the frequency 
of excess red-type Ha emitters relative to the total number 
of red-type stars. We restricted our analysis such that the 
magnitude of (B-V) colors used to identify candidate Be 
stars in each cluster (i.e. ~ -0.3 < (B-V) < 0.2 in NGC 
458, WB06) matched the range of (B-V) colors used to 
study the frequency of red-type excess emitters (i.e. 0.2 
< (B-V) < 0.7 for NGC 458). Furthermore, we adopted 
the same (R-Ha) cutoffs delineating excess emitters from 
normal stars for each cluster as used by WB06. The re- 
sults of this exercise are tabulated in column 2 of Table [T2l 
Interestingly, with few exceptions, the frequency of poten- 
tial "contaminants" implied by this technique is generally 
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consistent with the minimum rate of contamination de- 
rived from our polarization results (see column 3 of Table 
[T2l column 6 of Tabled, i.e. the candidates we assigned a 
Type 4 classification. Clearly this photometric contamina- 
tion check can not provide a diagnostic on the Be classifi- 
cation of individual objects, unlike our polarimetry; how- 
ever, we suggest that this technique might provide a rough 
estimate of the role of contaminants in cases in which de- 
tailed follow-up observations have not been made or are 
unfeasible. 

5.3. The Evolutionary Status of Excess Ha Emitters in 
Very Young (5-8 Myr) Clusters: Early-Type 
Classical Be stars or Early-Type Herbig Be stars? 

One of the interesting results presented in WB06 was 
the identification of candidate Be star in clusters of age 
5-8 Myr, in an abundance similar to the general frequency 
of the Be phenomenon in our Galaxy, ~17%. However, 
it was uncertain based on those results whether the de- 
tected excess Ha emitters a) were true classical Be star- 
disk systems; b) were B-type objects which still possessed 
remnant star-formation disks, i.e. Herbig Be stars; or c) 
merely had diffuse Ha nebulosity coincidentally associated 
with them. Our intrinsic polarimctric dataset affords us 
the opportunity to further constrain the evolutionary sta- 
tus of candidate Be stars in four clusters of age 5-8 Myr, 
LH 72, NGC 1858, NGC 346, and NGC 371. 

5.3.1. Intrinsic Polarization Signatures of Dust versus 

Gas Disks 

As discussed in the Introduction, Section [U and refer- 
ences therein, it is well accepted that classical Be stars 
are characterized by geometrically thin gas disks, and the 
mechanism responsible for producing their intrinsic polar- 
ization is electron scattering. Herbig Be stars are more 
complex systems in that their circumstellar environments 
contain both gas and dust. The dominant polarigenic 
mechanism in Herbig Be stars is well established to be scat- 
tering by dust, as out lined in the review papers oflBastienl 
(|l988f ) jGrinm1 (|l994f) and lWaters fe Waelkensl (|l998f) and 
numerous references therein. While Herbig Be and classi- 
cal Be stars exhibit some similar observational traits, such 
as the presence of an intrinsic, often variable, polariza- 
tion component (IVrba et al.lll979tlQuirrenbach et al.ll997t 
Vink et all l2002t iBiorkman fc Meadd 120051 : iBiorkmanl 



2006), the different scattering mechanisms which charac- 



terize Herbig Be versus classical Be stars manifest them- 
selves in several observational manners. The magni- 
tude of linear polarizati on of classical Be star s is typi- 
cally of order 1% or less (jBiorkman et al.l[20 00), whereas 
the average (~3%) and maximum (14.5%) polarization 
observed in Herbig Ae/Be stars is clearly much larger 
(jTamura fc Fu kagawa 2005). More importantly, the dif- 
ferent scattering mechanisms which dominate in each of 
these star-disk systems will produce a different wavelength 
dependence in these stars' intrinsic polarization. 

The characteristic wavelength dependent intrinsic po- 
larization signature of classical Be stars has already been 
documented in Section [4j Herbig Be stars are not known 
to exhibit any singular " characterist i c" intrinsic po l ariza- 
tion behavior (see e.g. IVrba et al.l [19791 iBastienl [1981 
iMever et all [2002); the grain chemistry, grain size distri- 
bution, and relative size and extent of the gas and dust 



disks/envelopes of individual Herbig Be systems likely con- 
tributes to the heterogeneous mix of observed intrinsic 
polarization. Some typical intrinsic polarization trends 
are observed in Herbig Ae/Be stars however, most no- 
tably rotations (and dramatic 90° reversals) in the in- 
trinsic polarization position angle (|Schmidt et alj [l992; 
IBiorkman et~all Il995l 119981 : IMever et alj l2002ft . We re- 
mark that several literature works which describe and in- 
terpret the wavelength dependence of intrinsic polariza- 
tion from Herbig Be stars, and at times attempt to draws 
parallels with the intrinsic polarization behavior of classi- 
cal Be stars, sho uld be v i ewed w ith ca ution. For e xam- 
ple, as noted by IBastienl (I1988D . both IVrbal (|l975l ) and 
iGarrison fc Andersonl (|l978l ) employ dubious assumptions 
in their efforts to constrain the interstellar polarization 
along the line of sight to their Herbig Be stars, and as such 
the wavelength dependence of the "intrinsic" polarization 
they report is likely inaccurate. Similarly, the discussion 
of the wavelength dependence of i ntrinsic polarizati on of 
the Herbig Be star BD+61 154 in IVrba et al.l (|1979D (see 
their Figure 2) ignores the wavelength dependence of the 
intrinsic polarization position angle. The ~40° intrinsic 
polarization position angle rotation from the B to I filter 
in these data is inconsistent with the expected signature of 
a pure gas disk system; the wavelength dependence of in- 
trinsic polarization from this Herbig Be star is not similar 
to that observed for classical Be stars. 

Although dust scattering is the dominant polarigenic 
agent for Herbig Be stars, it is still possible to detect evi- 
dence of the gaseous inner disk component of some of these 
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systems via polarimetry (Vink et all 120021 : 
2002| lOudmaiier et al.l l2005f) . Forexample. IMever et al. 
(|2002f ) presents intrinsic spectropolarimetry of the Herbig 
Be star MWC349A covering the wavelength range ~5000- 
10500 A (V,R,I filters), which exhibits both line depolariza- 
tion effects and a Paschen polarization jump (see their Fig- 
ure 7 ), similar to that observed in the classi ca l Be star zeta 
Tau (Wood. Biorkman. fc Biorkmanlll997l) . IMever et all 
(2002) suggest electron scattering as the mechanism re- 
sponsible for producing these specific features. Evidence 
of the dusty component of MWC349A's composite disk is 
also evident in these data. Specifically, the 30-40° rota- 
tion in the intrinsic continuum polarization position an- 
gle, beginning at ~6200 A and extending to the blue edge 
of the available data, ~5000 A indicates the presence of 
an additional scattering mechanism which becomes more 
pronounced at shorter wavelengths, i.e. dust scattering. 

Based on the available literature, we are aware of no 
compelling examples of Herbig Be stars whose intrinsic 
polarization magnitude and position angle fully follow the 
known UBVRI behavior of classical Be stars. While nei- 
ther typical nor likely, we can not rule out the possibility 
that one could in principle observe a Balmer jump sig- 
nature in the polarization magnitude of some Herbig Be 
stars, as these systems do have composite gas plus dust 
disks. However, as dust scattering is the dominant polar- 
igenic scattering agent in Herbig Be stars, it is dubious 
that the entire UBVRI intrinsic polarization magnitude 
and position angle of such systems would exhibit no evi- 
dence of a non-electron scattering component. Moreover, 
as astutely pointed out by our referee, it would be even 
more unlikely to observe one of our "Type 2 electron scat- 
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tering (ES) signatures" in any Herbig Be star, unless the 
dust grain chemistry and size distribution would conspire 
to produce grey scattering. 

Hence, the UBVRI spectral range of the intrinsic po- 
larization data presented in this paper enables us to dis- 
criminate deviations from a pure electron scattering signa- 
ture over a broad spectral range. This allows us to differ- 
entiate Herbig Be stars analogous to MWC 349A, whose 
gaseous disk dominates a portion of the optical polarimet- 
ric regime, and Herbig Be stars whose dust disks dominate 
the optical polarimetric regime, from classical Be stars. 

5.3.2. Nebular Origin? 

As noted in WB06, many of the candidate Be stars 
identified in clusters of age 5-8 Myr reside nearby or 
even within regions characterized by diffuse nebular emis- 
sion. We consider whether the observed characteristic po- 
larimetric signatures of classical Be stars, both electron 
scattering and Balmcr jump polarization features, could 
merely be artifacts of this background emission. Such po- 
larization signatures are only produced when the electron 
scattering optical depth is ~1.0, i.e., that present in the 
innermost disk region of classical Be stars. The density of 
any coincident diffuse Ha nebulosity present around can- 
didate Be stars is vastly insufficient to create such a po- 
larization signature. Thus while some of the "Be star" de- 
tectio ns made by Ha spectroscopic surveys ()Mazzali et al.l 
1996|) have been shown to be spurious detections due to the 
presence of coincident gas (|Keller fc Besselll 1 1 9981 ). such 
diffuse nebulosity can not be responsible for creating the 
polarimetric signals found in our dataset. 

5.3.3. Near-IR Colors 

Further insight into the evolutionary status of candi- 
date Be stars in our very young (5-8 Myr) clusters may 
be elicited via inspection of their near-IR colors, ft has 
been shown that the fundamentally different composition 
of classical Be disks (gas) and Herbig Ae/Be disks (gas 
and dust) will generally result in these objects occupying 
distinctly differen t locations in (J-H) versus (H-K) nea r-IR 
2-color diagrams (|Lada fc AdamJl992t iLi et al.lll99l , al- 
though minor instances of overlap are observed. As such, 
we cross-correlated all photometrically identified candi- 
date Be stars resi ding in clusters stud ied in this paper 
with the 2MASS (jSkrutskie et alj I2006D survey. Targets 
for which we were able to extract full JHK photometric 
data are listed in Table [TBI 

We plot the 2MASS JHK photometry of all available 
sources as a function of our polarimetric classifications in 
Figure [T5] (with error bars) and Figure [TBJ (without er- 
ror bars). Our Type-1 sources which exhibit polarimetric 
Balmer jumps are plotted as red circles, our Type-2 sources 
which exhibit ES polarimetric signatures are plotted as 
green triangles, and all other candidate Be stars (other 
Type-2 sources, Type-3 sources, Type-4 sources, and po- 
larimetric non-detections) are plotted as blue squares. For 
reference, we have also plotted the near -IR colors of 101 
known Galactic Be stars tabulated by iDoughertv et alj 
dl99l|) as filled bla c k tria ngles, 21 LMC "ELHC" stars 
which Ide Wit et a l. (2005) suggest are likely to be classi- 
cal Be stars as yellow open triang les, 2MASS and ground- 
based (jGummersbach et al.lll995l ) colors of 13 Magellanic 



Cloud B[e] stars known to have dust y disks as ligh t blue 
crosses, and the star "ELHC-7" which lde Wit et al.l (|2005h 
suggest is likely to be a LMC Herbig Ae/Be dusty star-disk 
system as a pink triangle. 

Figures [15] and [IB] clearly demonstrate that most of the 
candidate Be stars which we were able to correlate with 
the 2MASS catalog do lie in a distinctly different region 
of the JHK 2-color diagram than that occupied by known 
dustier (Herbig Ae/Be or B[e]) systems. Furthermore, al- 
though the Magellanic Cloud candidate Be stars' colors 
exhibit a wide dispersion, likely due in part to the consid- 
erable photometric errors of these faint sources (Figurc [T5l 
Table [TB]) . we find no clear distinction amongst the near- 
IR colors of our candidate Be stars as a function of their 
polarimetric classification. A small number of candidate 
Be stars marginally overlap with the location of dustier 
B[e] sources, although these sources are within 3-<r of the 
mean locus of classical Be colors. Clearly better quality 
near-IR photometry of these sources should be pursued 
to better elucidate their evolutionary status. Two can- 
didate Be stars clearly exhibit near-IR colors consistent 
with dusty disk systems, NGC 346:KWBBe 200 and NGC 
456WBBe 7. NGC 346:KWBBe 200 was not detected 
in our polarimetric survey owing to contam i nation from 
a nearby source; however, IWisniewski et al.l (|2007bl ) has 
recently published optical spectroscopic, near-IR photo- 
metric, and IR photometric observations of the target and 
concluded it is a B[e] supergiant, not a classical Be star. 
NGC 456 WBBe 7 was observed in our polarimetric survey 
and assigned a classification of Type-3. 5. Although noisy, 
our intrinsic polarization data exhibited marginal evidence 
of a position angle rotation, which could be indicative of a 
dusty envelope. The star's near-IR colors seems to verify 
this tentative interpretation. 

For comparison, we also plotted the same near-IR colors 
for only candidate Be stars residing in our very young (5-8 
Myr) clusters in Figure [171 Comparison of Figures [TBI and 
1 171 reveal no distinctive differences; rather, it is quite clear 
that the vast majority of candidate Be stars in these very 
young clusters have near-IR colors consistent with those 
expected from classical Be stars and notably inconsistent 
with those expected from composite gas plus dust disk 
systems (Herbig Ae/Be and B[e] stars). Similar to Figure 
[TBI Figure [T7] reveals that there is no clear color difference 
amongst our various polarimetric classifications. We in- 
terpret these results as additional evidence that the many 
of the candidate Be stars identified in clusters 5-8 Myr old 
by WB06 are classical Be stars, and not pre-main-sequence 
Herbig Be stars. 

5.3.4. Natal Disk Clearing Time 

The clearing time-scale of natal star formation material 
(gas and dust) will influence whether it is possible or likely, 
a-priori, to observe early-type Herbig star-disk systems in 
clusters of age 5-8 Myr. Amongst the less massive Herbig 
Ae and T Tauri stars, the median lifetime of inner optically 
thick accretion disks is 2-3 Myr, and although exceptions 
do exist, there is little to no H- and K-band evidence of 
primordial disks b eyond a median stellar age of 5 Myr 
(Hillenbrand 2005). The typical natal disk dissipation for 
more massive stars is less well known. From a detailed 
study of the young cluster NGC 6611, Ide Winter et al.l 



11 



(|1997| ) suggested that the clearing of natal disk material 
around the more massive st ars in this cluster was typically 
0.1 Myr or less. Similarly. iPogodin et a l. (2006) recently 
suggested that BOe star HD 53367 is a classical Be star, 
based in part on an observed episodic loss of its disk mate- 
rial, a ph enomenon known to cha ra cterize many c l assica l 
Be stars (<Porter fc Riviniusll2003l ). IPogodin et all (|2006l) 
also identified a 4-5 M so i ar pre-main-sequence binary com- 
panion which, given the pre-m ain-sequence evolutionary 
tracks of iPalla fc Stahlerl (|l993), suggests the system's age 
is < 0.8 Myr and indicates the natal gas plus dust disk of 
the BOe star HD 53367 must have dissipated at least on 
this time-scale. Although clearly not conclusive, these re- 
sults do suggest that it is reasonable to expect that most 
of the natal disks will have been cleared from the early 
B-type stars in our clusters of age 5-8 Myr. 

5.3.5. Summary: Most Candidate Be stars in Clusters 
5-8 Myr Old Are Classical Be Stars 

WB06 and iKeller et"ail (|l999t) identified a large pop- 
ulation of B-type stars in clusters of age 5-8 Myr which 
exhibited an excess of Ha emission. Observational studies 
suggest that the natal star formation disks of early B-type 
stars can clear on timescales of at least 0.1-0.8 Myr, al- 
though it is admittedly not clear if this timescale is "typi- 
cal" for all early-type B stars. The near-IR colors of many 
of the excess Ha emitters in these very young clusters arc 
consistent with those expected from gaseous classical Be 
stars and generally inconsistent with those observed from 
dustier B[e] and/or Herbig Ae/Be disk systems. The ob- 
served wavelength dependence of intrinsic polarization of 
many Ha excess stars in clusters of age 5-8 Myr are also 
consistent with that expected from pure gas disk systems 
(classical Be stars) and exhibit no evidence of secondary 
contributions from dust scattering, as would be expected 
for composite disk systems (Herbig Ae/Be, B[e] stars). 
Based on these factors, we suggest that there is compelling 
evidence supporting the existence of classical Be stars in 
clusters of age 5-8 Myr. 



5.3.6. Be Stars in the First Half of a B Star's Main 

Sequence Lifeti 

VFabreaat & Torreion 



Sequence Lifetime: A Co mparison to 

200(1) and \Martavan et al\ 



As outlined and summarized in Section [5. 3. 5[ data pre- 
sented in this paper support the suggestion of WB06 that 
classical Be stars do develop before the mid-point main 
sequence lifetime of a B star. Specifically, these data sup- 
port the presence of a significant number of classical Be 
stars, having crude photometric spectral types of B0-B5 
(WB06), in clusters spanning a near-continuous range of 
ages from log(t) of 6.7-8.1, i. e. 5-126 Myr. These results 
conf lict with claims made b y lFabregat fc Torreion! (2000) 
and iMartavan et alj {2007) that the Be phenomenon de- 
velops in the second-half of a B star's main sequence 
lifetime. IMartavan et alj (|2007h do suggest that mas- 
sive Be stars may briefly appear near the ZAMS in the 
SMC and LMC, before losing their "Be status" until 
the very end of the first part of their main sequence 
lifetimes. We note that this conclusion was essentially 
based on a single cluster population, which limits its ro- 
bustness. WB06 identified, and the intrinsic polariza- 
tion presented in this paper support, the detection of a 



large body of Magellanic Cloud Be stars having crude, 
photometrically-assigned spectral types of B0-B5, in clus- 
ters spanning a range of ages from 5-126 Myr. No abrupt 
absence of classical Be stars was observed in any of the 
younger- type clusters of our photometric and/or polari- 
metric datasets. Thus, our larger, more complete inves- 
tigation of Be stars in Magellanic Cloud clusters calls 
into question the suggestion that the Be phenomenon 
is primarily restricted to th e second-half of the main- 
sequence lifetime of B -stars (jFabregat fc Torreion! 120001 : 
IMartavan et a l. 2007) and contradicts the appearance- 
disappear ance-reappearance of t he Be phenomenon sug- 
gested by IMartavan et all (|2007t) . 

5.3.7. Implications of the Presence of Be Stars in Very 
Young Clusters 

WB06 first reported the presence of candidate classical 
Be stars in clusters as young as 5 Myr, and noted that if 
such objects were truly classical Be stars, they would not 
have spent enough time on the main sequence to spin-up 
to near-c ritical rotation velocities via the mechanism pro- 
posed bv lMevnet fc Maederl (I2000T) or via mass transfer in 
a binary system ( McSwain fc G ics 2005). The confirma- 
tion of the classical Be status of many of these stars via the 
present study can be interpreted as evidence that a signif- 
icant number of classical Be stars emerge onto the zero- 
age-main-sequence at near critical rotation velocities. An 
alternate interpretation of our results is that these youthful 
classical Be stars are evidence of the existence of a subset 
of the Be phenomenon which rotate at significantly sub- 
critical r ates, perhaps a s l ow as 0.4-0.6 \ C r it for early- type 
B stars (|Cranmerll2005[) . ICranmerl (|2005h suggested that 
only a subset of early-type (07-B2) Be stars might be sub- 
critical rotators, while later-type objects (B3-A0) should 
be all near-critical rotators. Our photometric survey of 
these very young clusters assigned crude spectral types to 
these candidate Be stars of B0-B5 (WB06); however, ow- 
ing to detection biases in our follow-up polarimetric ob- 
servations, we have in general only been able to confirm 
that portions of the brightest of these candidates (i.e. the 
earliest spectral sub-types, B0 to ^B3) are bona-fide clas- 
sical Be stars. Thus our present observations do not rule 
out the possibility that some of our very young bona-fide 
classical Be stars might belong to t he su b-critical rotation 
population predicted by ICranmerl (|2005l ) . 

The observed emergence of the Be phenomenon earlier 
in the main sequence lifetime than previously thought also 
has important implications regardi ng the role of magneti c 
fields in the formation o f Be d isks (|CassineHi et al.ll2003l ). 
iMacGregor fc Cassinellil ([2003) investigated the transport 
of magnetic flux tubes in 9 M su „ stars (^B2.5-B3) and 
found these structures could rise from the core to the sur- 
face within the (at the time) expected beginning of the Be 
phenomenon, the mid-point main sequence lifetime. The 
present study has confirmed the presence of bona-fide clas- 
sical Be stars in clusters as young as 5 Myr having crude, 
photometrically assigned spectral types of B0 to ^B3. Al- 
though WB06 suggested the presence of an additional pop- 
ulation of B4-B5 type candidate Be stars in these clus- 
ters, the limited dynamical range probed by the present 
study was unable to confirm the status of these fainter 
candidates. Although the magnetic flux transport model 
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of iMacGregor fc Cassinellil (|2003l ) had sufficient time to 
transport flux to the stellar surface for ^B3 stars, this is 
unlikely to be the case for later type stars (>B4) which 
have a much thicker radiative envelope. Hence, the ex- 
istence of such later-type Be stars in extremely young 
clusters would likely require an additional mechanism to 
be employed to accelerate the rise times of magnetic flux 
(Cassinelli, 2006 personal communication). While chal- 
lenging, follow-up spectroscopy of these young cluster pop- 
ulations would p rovide one avenue confirm the presence of 
these purported (jWisniewski fc Biorkmanl2006l ) later-type 
systems. 

5.4. Future Work 

While the present study provides significant advances 
in identifying and understanding the biases present in ear- 
lier 2-CD studies of cluster populations, additional work 
is clearly required. Our polarimctric survey was not sensi- 
tive to pole-on or near-pole-on systems, hence determining 
which of our "type-2" objects are near-pole-on classical Be 
stars and which are astrophysical objects of a fundamen- 
tally different nature is a high priority. We suggest that 
followup moderate resolution optical spectroscopic or in- 
frared photometric observations would provide reasonable 
diagnostics to resolve this bias, and enable more quanti- 
tative determinations of the bona-fide classical Be content 
of clusters as a function of age and/or metallicity. Our 
data suggest that at least 25% of photometrically identi- 
fied candidate Be stars appear unlikely to be true classical 
Be stars; however, we are unable to place firm constraints 
on the true astrophysical nature of such objects. We sug- 
gest that followup infrared photometric and/or spectro- 
scopic observations would be useful to further constrain 
the evolutionary nature of these objects. Furthermore, 
such observations would also be useful to ascertain the 
true nature of candidate Be stars which were too faint to 
be reliably probed by our polarimctric observations. 

Our analysis of the intrinsic polarization properties of 
LMC and SMC classical Be stars suggests that the funda- 
mental disk properties of classical Be stars may depend on 
metallicity. Modeling followup moderate resolution spec- 
tropolarimetric observations of LMC/SMC classical Be 
stars is a clear avenue one could use to further investigate 
these results, and we note that the advent of accurate lin- 
ear spectropolarimctcrs on large aperture telescopes (i.e. 
the Focal Reducer Spectrograph (FORS) at the Very Large 
Telescope (VLT) or the Robert Stobie Spectrograph at the 
South African Large Telescope (SALT)) make obtaining 
such observations feasible. 

6. SUMMARY 

We have used imaging polarimetric observations of six 
LMC and six SMC clusters to investigate the evolutionary 
status of B-type stars identified as excess Ha emitters via 
optical 2-color diagram photometric techniques. We char- 
acterized the interstellar polarization components of these 
data in a systematic manner, allowing us to isolate the in- 
trinsic polarization properties of our data, hence constrain 
the dominant polarigenic agent in each system. We found: 

1. The interstellar polarization associated with NGC 
330 was characterized by \ ma x ~4500A, suggesting 
the presence of a small dust grain population; 



2. The ISP of NGC 2100 exhibited clear evidence of 
a complex morphology, indicating the presence of 
a non-uniform magnetic field. We offer a detailed 
discussion of this B field and its potential origin in 
IWisniewski et all (|2007al ). 

3. The UBVRI wavelength dependence of intrinsic po- 
larization of many candidate Be stars in clusters 
of age 5-8 Myr exhibit either polarization Balmer 
jumps and/or electron scattering signatures, which 
is the expected diagnostic of free-free scattering 
from a gaseous disk. No evidence of a secondary 
polarigenic agent, i.e. dust scattering, is observed 
in these system. Moreover, the 2MASS near-IR col- 
ors of many of these systems are most consistent 
with the expected colors of pure gas disks, and in- 
consistent with the observed colors of dustier Mag- 
ellanic Cloud disk systems (Herbig Ae/Be, B[e]). 
We con clude these data confirm t he ini tial sugges- 
tion of IWisniewski fc Biorkmanl (|2006[ ) that clas- 
sical Be stars are present in clusters of age 5-8 
Myr, contradicting claims that the Be phenomenon 
only develops in the second half of a B star's 
mam sequence l ifetim e (|Fabregat fc Torre ionll2000l; 
iMartavan et~aT]|2007l) . i.e. after 10 Myr. 

4. We interpret the observed presence of classical Be 
stars in clusters of age 5-8 Myr as evidence that 
a significant population of early-type B stars must 
emerge onto the ZAMS rotating at near critical ve- 
locities. We note however that these results do not 
exclude the possibility that we are observing the hy- 
pothesized subset of clas sical Be stars w hich rotate 
at sub-critical velocities ()Cranme n l2005h . 

5. Comparing the polarimetric properties of our 
dataset to a similar survey of Galactic classical Be 
stars, we find the prevalence of polarimetric Balmer 
jump signatures decreases with metallicity. We 
speculate that these results might indicate that ei- 
ther it is more difficult to form large disk systems 
in low metallicity environments, or that the average 
disk temperature is higher in these low metallicity 
environments. 

6. We find evidence that at least 25% of photometri- 
cally identified candidate Be stars do not exhibit po- 
larimetric signatures consistent with those expected 
from classical Be stars. These data strongly suggest 
that caution must be exercised when attempting 
to correlate the onset of the Be phenomenon with 
evolutionary age and/or metallicity based solely on 
statistics derived from simple 2-color diagram pho- 
tometry. 
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Table 1 
Summary of our observations 



Cluster 


Location 


Filter 


Age 


Date 




Exposure Time 


Brack 60 


SMC 


u 


o 


2002 


Oct 


19 


1200 






B 




2002 


Oct 


21 


720 






V 




2002 


Oct 


19 


300 






R 




2002 


Oct 


19 


300 






I 




2002 


Oct 


21 


300 


NGC 330 


SMC 


U 


y 


2001 


Nov 


22 


1200 






B 




2001 


Nov 


23 


420 






V 




2001 


Nov 


25 


300 






R 




2001 


Nov 


25 


240 






I 




2001 


Nov 


24 


240 


NGC346 


SMC 


U 


vy 


2001 


Nov 


21 


720 






B 


2001 


Nov 


22 


240 






V 




2001 


Nov 


21 


180 






R 




2001 


Nov 


24 


180 






I 




2001 


Nov 


23 


180 


NGC 371 


SMC 


U 


vy 


2002 


Oct 


18 


1200 






B 




2002 


Oct 


18 


600 






V 




2002 


Oct 


17 


420 






R 




2002 


Oct 


17 


300 






I 




2002 


Oct 


18 


420 


NGC 456 


SMC 


U 


y 


2002 


Oct 


24 


1200 






B 




2002 


Oct 


26 


720 






V 




2002 


Oct 


24 


420 






R 




2002 


Oct 


26 


300 






I 




2002 


Oct 


25 


300 


NGC 458 


SMC 


U 


o 


2002 


Oct 


27 


1200 






B 




2002 


Oct 


27 


720 






V 




2002 


Oct 


26 


300 






I 




2002 


Oct 


27 


300 


LH 72 


LMC 


u 


vy 


2002 


Oct 


26 


1200 






B 




2002 


Oct 


25 


720 






V 




2002 


Oct 


25 


420 






R 




2002 


Oct 


27 


360 






I 




2002 


Oct 


26 


300 


NGC 1818 


LMC 


U 


y 


2001 


Nov 


21 


900 






B 




2001 


Nov 


22 


240 






V 




2001 


Nov 


21 


240 






R 




2001 


Nov 


22 


180 






I 




2001 


Nov 


22 


300 


NGC 1858 


LMC 


U 


vy 


2002 


Oct 


19 


1200 






B 




2002 


Oct 


24 


720 






V 




2002 


Oct 


18 


420 






R 




2002 


Oct 


21 


300 






I 




2002 


Oct 


24 


300 


NGC 1948 


LMC 


U 


y 


2001 


Nov 


27 


1200 






B 




2001 


Nov 


26 


360 






V 




2001 


Nov 


25 


180 






R 




2001 


Nov 


26 


240 






I 




2001 


Nov 


26 


300 


NGC 2004 


LMC 


U 


y 


2001 


Nov 


25 


900 






B 




2001 


Nov 


22 


180 






V 




2001 


Nov 


27 


180 






I 




2001 


Nov 


22 


420 


NGC 2100 


LMC 


u 


y 


2001 


Nov 


24 


1200 






B 




2001 


Nov 


23 


240 






V 




2001 


Nov 


23 


180 






R 




2001 


Nov 


24 


180 






I 




2001 


Nov 


23 


180 



Note. — We summarize some of the properties of the data presented in this paper. The exposure times lis ted are the integration time use d 
at each of eight wave-plate position. The age labels in column 4, originally defined and discussed in detail in Wisnicwski & Biorkman (2006), 
correspond to very young (vy) 5-8 Myr; young (y) 10-25 Myr; and old (o) 32-158 Myr. 
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Table 2 

Average interstellar polarization 



Cluster 


% P„ 


% Pt, 


% P„ 


% P r 


% P, 


PA 


P m „ 




dPA 


K 


Bruck 60 


0.46±0.187o 


0.51±0.30% 


0.55±0.23% 


0.54±0.19% 


0.48±0.19% 


130 


0.55% 


5500 





0.923 


LH 72 


0.34±0.14% 


0.37±0.12% 


0.40±0.14% 


0.39±0.12% 


0.35±0.14% 


28 


0.40% 


5500 





0.923 


NGC 330 


0.51±0.17% 


0.53±0.16% 


0.52±0.14% 


0.48±0.17% 


0.41±0.17% 


120 


0.53% 


4500 





0.737 


NGC 346 


0.31±0.17% 


0.35±0.18% 


0.37±0.17% 


0.36±0.15% 


0.32±0.19% 


125 


0.37% 


5500 





0.923 


NGC 371 


0.38±0.15% 


0.42±0.16% 


0.45±0.15% 


0.44±0.15% 


0.39±0.15% 


118 


0.45% 


5500 





0.923 


NGC 456 


0.57±0.21% 


0.63±0.24% 


0.67±0.25% 


0.66±0.29% 


0.58% 


147 


0.67% 


5500 


-8 


0.923 


NGC 458 


0.33±0.34% 


0.36±0.23% 


0.39±0.22% 




0.34±0.15% 


120 


0.39% 


5500 





0.923 


NGC 1818 


0.52±0.17% 


0.58±0.16% 


0.62±0.15% 


0.61±0.17% 


0.54±0.14% 


39 


0.62% 


5500 





0.923 


NGC 1858 


0.32±0.21% 


0.36±0.15% 


0.38±0.16% 


0.37±0.22% 


0.33±0.17% 


45 


0.38% 


5500 





0.923 


NGC 1948 


0.57±0.19% 


0.64±0.15% 


0.68±0.18% 


0.67±0.17% 


0.59±0.17% 


35 


0.68% 


5500 





0.923 


NGC 2004 


0.32±0.12% 


0.36±0.13% 


0.38±0.15% 




0.33±0.16% 


26 


0.38% 


5500 





0.923 



Note. — The Serkowski parameters PA, Pmax, ^max, dPA, and K which best described each of our clusters are tabulated, along with the 
final ISP values in the U, B, V, R, and I filters. The horizontal line separates SMC clusters (above) from LMC clusters (below). 



Table 3 
Foreground ISP Corrections 



Schmidt Region Polarization Position Angle Clusters 

Region I 0.37 ± 0.15% m Bruck 60 

Region II 0.27±0.15% 123 NGC 330, 346, 371 

Region III 0.19±0.21% 110 NGC 456, 458 

Region VI 0.64±0.19% 30 LH 72 

Region VII 0.32±0.13% 29 NGC 1818,1858,1948 

Region IX 0.40±0.13% 20 NGC 2004 



Note. — We used the Galactic polarization maps of[Schmidt (1976|) as templates to identify and remove the foreground Galactic interstellar 
polarization contribution to each of our Magellanic Cloud lines of sight. Columns two and three represent the corrections applied to each of 
our lines of sight, listed in column four. 



Table 4 

ISP Intrinsic to the Magellanic Clouds 



Cluster 


P u 


PA„ 


Pi, 


PA b 


Pi, 


PA„ 


Pr 


PA r 


Pi 


PA, 


Bruck 60 


0.43±0.23% 


132±15 


0.45±0.34% 


133±22 


0.34±0.27% 


151±23 


0.41±0.24% 


138±17 


0.33±0.24% 


144±21 


LH 72 


0.27±0.24% 


28±26 


0.24±0.22% 


27±26 


0.24±0.24% 


123±29 


0.07±0.22% 


18±90 


0.07±0.24% 


130±98 


NGC 330 


0.48±0.23% 


126±14 


0.48±0.22% 


126±13 


0.25±0.21 


129±24 


0.34±0.23% 


127±19 


0.24±0.23% 


128±28 


NGC 346 


0.28±0.23% 


125±24 


0.30±0.23% 


125±22 


0.10±0.23% 


130±66 


0.22±0.21% 


126±27 


0.15±0.24% 


127±46 


NGC 371 


0.35±0.21% 


118±17 


0.37±0.22% 


117±17 


0.19±0.21% 


111±32 


0.31±0.21% 


116±19 


0.22±0.21% 


114±27 


NGC 456 


0.56±0.30% 


148±15 


0.62±0.32% 


149±15 


0.64±0.33% 


155±15 


0.64±0.36% 


151±16 


0.56% 


153.2 


NGC 458 


0.31±0.40% 


121±37 


0.32±0.31% 


121±28 


0.22±0.30% 


129±39 


0.23±0.26% 


125±32 






NGC 1818 


0.49±0.21% 


40±12 


0.52±0.21% 


40±12 


0.34±0.20% 


49±17 


0.46±0.21% 


43±13 


0.35±0.19% 


45±16 


NGC 1858 


0.29±0.25% 


47±25 


0.31±0.20% 


48±19 


0.20±0.21% 


74±30 


0.25±0.26% 


55±30 


0.19±0.21% 


63±32 


NGC 1948 


0.54±0.23% 


35±12 


0.58±0.20% 


36±10 


0.37±0.22% 


40±17 


0.51±0.21% 


37±12 


0.39±0.21% 


38±15 


NGC 2004 


0.28±0.18% 


27±18 


0.28±0.18% 


28±18 


0.08±0.20% 


75±72 


0.09±0.21% 


44±67 







Note. — We corrected the total measured ISP components of each of our clusters for the effects of foreground Galactic interstellar polarization, 
revealing the ISP components which arose from the dichroic absorption of light by dust grains located within the Magellanic Clouds. The 
horizontal line separates SMC clusters (above) from LMC clusters (below). 
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Table 5 

Total interstellar polarization toward NGC 2100 



Cluster 




% P„ 


% P„ 


% P„ 


% P r 


% P* 


PA 


P m ax 




dPA 


K 


NGC 2100 


area 1 


1.32% 


1.41% 


1.44% 


1.37% 


1.20% 


139 


1.45% 


5000 





0.83 


NGC 2100 


area 2 


1.26% 


1.35% 


1.37% 


1.31% 


1.14% 


100 


1.38% 


5000 





0.83 


NGC 2100 


area 3 


1.64% 


1.69% 


1.65% 


1.54% 


1.32% 


78 


1.70% 


4500 





0.737 



Note. — A summary of the total interstellar polarization values toward NGC 2100 for the 3 distinct spatial regions identified within the 
field of view of our NGC 2100 dataset. 



Table 6 

Total polarization of candidate Be stars 



Candidate Be Star 


Filter 


%p 


PA 


% Q 


% u 


% Err 


Bruck 60:WBBe 1 


u 
b 


1.75 


107.3 


-1.44 


-0.99 


0.61 




V 


1.33 


113.7 


-0.90 


-0.98 


0.30 




r 


1.33 


117.6 


-0.76 


-1.09 


0.35 



Note. — The total polarization (i.e. interstellar plus intrinsic components) of one candidate Be star investigated in this study is tabulated. 
A table presenting results for our entire dataset is available in the electronic version of this Journal. 



Table 7 

Intrinsic polarization of candidate Be stars 



Candidate Be Star 


Filter 


%p 


PA 


% Q 


% U 


% Err 


Bruck 60:WBBc 1 


u 
b 


1.46 


100.9 


-1.36 


-0.54 


0.61 




V 


0.91 


104.4 


-0.80 


-0.44 


0.30 




r 


0.88 


110.0 


-0.67 


-0.56 


0.35 



Note. — The intrinsic polarization of one candidate Be star investigated in this study is tabulated. A table presenting results for our entire 
dataset is available in the electronic version of this Journal. 



Table 8 

Classifying the intrinsic polarization of individual candidate Be stars 

Name Pol Balmcr J umpY PA Classification Comments 

Bruck 60:WBBc 1 N PI 2 moderate ES 



Note. — The classification of the intrinsic polarization signature of one of the candidate Be stars investigated in this study is presented. 
Column 2 indicates whether a polarization Balmcr jump was present in the data; the letters in column 3 indicate whether the polarization 
position angle was fiat (F), linearly sloped (S), or randomly variable (V) while the numbers indicate the a level of such characterizations; 
column 4 indicates the likelihood each object is a bona-fidc classical Be star based on the 4 point classification scheme described in the text; 
and any relevant comments are listed in column 5. Note that a table presenting results for our entire dataset is available in the electronic 
version of this Journal. 
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Table 9 

Intrinsic polarization summary 



Cluster 


Octcction Rate 


# Pol. BJ 

(type-1) 


# ES 
(type-2) 


Not Inconsistent 
(type-2) 


Unlikely 
(type-3,4) 


Bruck 60 
NGC 330 
NGC 346 
NGC 371 
NGC 456 
NGC 458 


18/26 (69%) 
41/76 (54%) 
33/48 (69%) 
73/129 (57%) 
14/23 (61%) 
10/30 (33%) 


1/18 (6%) 
0/41 (0%) 
8/33 (24%) 
10/73 (14%) 
0/14 (0%) 
0/10 (0%) 


8/18 (44%) 
9/41 (22%) 
7/33 (21%) 
6/73 (8%) 
1/14 (7%) 
0/10 (0%) 


11/18 (61%) 
24/41 (59%) 
22/33 (67%) 
49/73 (67%) 
9/14 (64%) 
6/10 (60%) 


6/18 (33%) 
17/41 (41%) 

3/33 (9%) 
14/73 (19%) 
5/14 (36%) 
4/10 (40%) 


LH 72 
NGC 1818 
NGC 1858 
NGC 1948 
NGC 2004 
NGC 2100 


34/50 (67%) 
18/40 (45%) 
27/39 (69%) 
22/27 (81%) 
43/67 (64%) 
35/61 (57%) 


1/34 (3%) 
4/18 (22%) 
3/27 (11%) 
6/22 (27%) 
9/43 (21%) 
8/35 (23%) 


5/34 (15%) 
5/18 (28%) 
2/27 (7%) 
6/22 (27%) 
10/43 (23%) 
16/35 (46%) 


22/34 (65%) 
12/18 (67%) 
13/27 (48%) 
12/22 (55%) 
28/43 (65%) 
20/35 (57%) 


11/34 (32%) 
2/18 (11%) 
11/27 (41%) 
4/22 (18%) 
6/43 (14%) 
7/35 (20%) 



Note. — We summarize the statistics of our classification of the intrinsic polarization properties of our dataset. The detection rate column 
describes the number of candidates for which we were able to extract polarimetric information. The third column denotes the fraction of 
candidates which display polarimetric Balmer jumps (BJ): we consider these objects to be bona-fide classical Be stars. The fourth column 
describes the fraction of candidates which exhibit electron scattering (ES) signatures, i.e. moderate amounts of wavelength-independent 
polarization at a wavelength independent position angle. The fifth column describes candidates that, to within 3 <r, can not be ruled out as 
possible Be stars. The sixth column describes candidate Be stars whose intrinsic polarization docs not agree, to within 3 <r, with that expected 
from classical Be stars. The horizontal line separates SMC clusters (above) from LMC clusters (below). 



Table 10 

Average intrinsic polarization properties 



Cluster 


Mean BJ 


Median BJ 


Mean ES 


Median ES 


Mean Unlikely 


Median Unlikely 


Very Young 
Young 


13% 
13% 


13% 
21% 


13% 
22% 


12% 

23% 


25% 
26% 


26% 
20% 


SMC 
LMC 


7% 
18% 


3% 
22% 


17% 

24% 


15% 

25% 


30% 
23% 


35% 
19% 


ALL SMC+LMC 


13% 


13% 


21% 


22% 


26% 


26% 



Note. — Statistical averages of data listed in Table [9] are presented. Clusters were grouped according to age (very young and young) as well 
as metallicity (SMC and LMC). The final entry labeled ALL denotes a global average of all SMC and LMC cluster data regardless of age or 
metallicity. 



Table 11 

Average intrinsic polarization properties with "contaminants" removed 



Cluster 


Mean BJ 


Median BJ 


Mean ES 


Median ES 


Very Young 


17% 


18% 


17% 


18% 


Young 


16% 


24% 


28% 


31% 


SMC 


7% 


4% 


25% 


17% 


LMC 


22% 


25% 


31% 


29% 


ALL LMC+SMC 


16% 


18% 


28% 


25% 


Milky Way Galaxy 


42-49% 









Note. — Objects deemed unlikely to be classical Be stars in Table [TOl were removed from consideration, and the statistical averages of data 
listed in Table|9]were recalculated. The "ALL SMC+LMC" entry denotes a global average of all SMC and LMC cluster data regardless of age 
or metallicity. The Milky Way data averages were derived from analysis of the HPOL polarization catalog, as discussed in Section 5.1. 
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Table 12 

Frequency of Contaminants in Be Populations Identified via 2-CDs 

Cluster Analysis oi 2-CD data Analysis oi Polarimctry 

Brack 60 14% 33% 

NGC 371 11% 19% 

NGC 456 10% 36% 

NGC 458 38% 40% 

LH 72 25% 32% 

NGC 1858 42% 41% 



Note. — The possible rate of contamination of classical Be star detections extracted from photometric 2-color diagrams. As discussed in 
Section 5.2, the second column represents the frequency of red-type excess Ha emitters, which we suggest might be an appropriate proxy of 
the frequency of B-type excess Ha emitters which are not classical Be stars. The third column is a reproduction of column 6 in Table [9] and 
illustrates the frequency of 2-CD contaminants as derived from polarimetric observations of candidate Be stars. 



Table 13 
IR Colors of Candidate Be Stars 



Star 


Classification 


(J-H) 


(H-K) 


(J-H) err 


(H-K)„ r 


Bruck60:WBBc 5 


Typc2-4, other 


0.09 


0.04 


0.14 


0.24 


Bruck60:WBBc 20 


Typc2-4, other 


0.29 


0.00 


0.13 


0.20 


LH72:WBBc 5 


Typel:Bal-J 


0.07 


0.21 


0.10 


0.15 


LH72:WBBc 15 


Typc2:ES 


0.54 


0.24 


0.14 


0.17 



Note. — Near-IR colors are tabulated for all candidate Be stars in clusters investigated in this paper which were detected in all three 
2MASS photometric bands. The quoted color errors were obtained by adding individual filter errors in quadrature. Column 2 designates the 
polarimetric designation, if any, assigned to each target within this paper; "Type-l:Bal-J" designates Type-1 sources which we assert are most 
likely to be classical Be stars, "Type-2:ES" designates Type-2 sources which exhibit electron scattering intrinsic polarimetric signatures and are 
thus likely classical Be stars, and "Type-2-4, other" designates all other sources, regardless of whether they were detected via our polarimetric 
survey or not. Note that a table presenting results for our entire dataset is available in the electronic version of this Journal. 
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Fig. 1. — The output of the PCCDPACK routine select for the V-band observation of NGC 1818 is shown. Cluster members having a 
polarimetric signal-to-noise level, p/cr p , greater than 5 are plotted on a Stokes Q-U diagram (panel a), as a polarization map as a function of 
pixel location (panel b), and in position angle (panel c) and polarization (panel d) histograms. As many cluster stars in panels a, c, and d 
exhibit the same general polarimetric properties, we suggest that most of these objects lack significant intrinsic polarization components and 
instead only exhibit evidence of interstellar polarization. 
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Fig. 2. — The observationally derived tot al ISP estimates for ea ch filter of NGC 1818 are plotted as a function of wavelength. The data 
closely follow a classic Serkowski-law shape llSerkowski et ahl[l"9 75), sugge sting that our method for determining the ISP in each cluster was 
reasonable. Over-plotted on these estimates is a modified Serkowski-law (Scrkowski ct al. 1975: Wilkine ct al. 1982) which wc deemed best 
fit the estimates. We extracted the U, B, V, R, and I filter polarization associated with this fit, and hereafter use these values to describe the 
ISP of NGC 1818. 




Fig. 3. — The observationally derived total IS P estimates for each fi l ter of NGC 330 are plotted as a function of wavelength. Over-plotted 
on these estimates is a modified Serkowski-law (Scrkowski ct al. 1975; Wilking ct al. 1982) which wc deemed best fit the data. These data, 
as well as the SMC ISP associated with NGC 330 (see Table |4j are characterized by a X m ax value of ~4500 A, suggesting the presence of 
smaller than average dust grains. 
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Fig. 4. — The I-band total polarization vectors of our observations of the LMC cluster NGC 2100 arc over-plotted on a DSS-2 red image. 
The general morphology of these polarization vectors matches that seen in the other 4 filters observed, indicating this structure is real. 



.Q JO- 
'S 

I 



NGC 2100 B 5% 

~ i — i — i — I — i — i — i — i — I — i — i — I — i — I — i I- .j| 



'".V" 



i i i i I i i i i I i i i i L_ 



5 b 4? 



4f | 

Right Ascension 



Fig. 5. — The B-band polarization vector map of region 1 in NGC 2100 ISP space is over-plotted on a DSS-2 blue image. 
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Fig. 6. — The B-band polarization vector map of region 2 in NGC 2100 ISP space is over-plotted on a DSS-2 blue image. 
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Fig. 7.— 



The B-band polarization vector map of region 3 in NGC 2100 ISP space is over-plotted on a DSS-2 blue image. 



23 



o 




175.0 . 



M 150.0 . 

0) 

Q 



PU 



125.0 



100.0 I 



4000, 



6000, 8000 

WAVELENGTH (A) 



10000, 



Fig. 8. — The observationally derived initial estimates of the total interstellar polarization towards region 1 of NGC 2100 is shown, along 
with the modified Serkowski-law we deemed best fit the data. The final ISP values for region 1 used in the remainder of this paper were 
extracted from this Scrkowski curve. 
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Fig. 9. — The time variability of the intrinsic polarization of the known Galactic classical Be star pi Aquarii is shown in these multi-epoch 
data from the HPOL spectropolarimcter. During the 1989 and 1991 observations, the inner disk density was sufficiently high such that the 
effects of pre- and post-scattering absorption by hydrogen in the disk is seen in the appearance of a polarization Balmer jump and a "saw- 
tooth" like wavelength dependence. The inner disk density was sufficiently low in the 1990 observation such that the wavelength independent 
electron scattering signature was not noticeably altered by a hydrogen opacity signature. We have examined our intrinsic imaging polarization 
observations of candidate Be stars to search for the polarimetric signatures found in this figure. 
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Fig. 10. — The intrinsic polarization of NGC 371:WBBcl8 exhibits the "sawtooth-like" polarization signature characteristic of classical Be 
stars. Such objects were assigned a designation of type-1 to indicate they are definitely classical Be stars. 
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Fig. 11. — The intrinsic polarization of NGC 371:WBBe21 clearly exhibits a wavelength independent electron scattering polarization 
signature. Such objects were assigned a designation of type-2 to indicate they are not inconsistent with being classical Be stars. 
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Fig. 12. — The intrinsic polarization of NGC 371:WBBe47 is not consistent, to within 3 cr, with that expected from a classical Be 
star. Such objects would receive a designation of typc-3 or typc-4, depending on the severity of their deviation from any of the expected 
polarimctric signatures. These objects (especially those which receive a type-4 designation) provide strong evidence that the 2-CD technique 
can mis-identify classical Be star-disk systems. 



27 



Q_ 

s 




< 
Q_ 



3000 5000 7000 9000 
Wavelength 



60 



40 



20 








3000 5000 7000 9000 
Wavelength 



Fig. 13. — The intrinsic polarization of NGC 1948:KWBBe 246 follows a complex wavelength dependence. The optical polarization 
appears to decrease and pos ition angle increase at short wavelengths. Objects which have dusty plus gaseous disks, such as HD 45677 
( Schultc-Ladbcck ct al. 1992), exhibit similar optical polarimetric behavior, suggesting that NGC 1948:KWBBe 246 and similar objects in 
our samp le might also be chara cterized by such disks. Alternatively, the wavelength dependence of these data also resemble a Serkowski-likc 
behavior (Scrkowski ct al. 1975), which as discussed in Section 4 might indicate that some these types of objects are characterized by abnormal 
interstellar dust conditions. 
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Fig. 14. — The wavelength dependent intrinsic polarization of NGC 2100:KWBBe 111 clearly is not consistent with that expected from a 
classical Be star. We suggest that this object might be a dust-disk system. 
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Fig. 15. — This near-IR 2-color diagram presents all available 2-MASS photometry of candidate Be stars residing in clusters investigate 
in this paper. "Type-l:Bal-J" sources (red circles) correspond to our polarimetrically designated Type-1 sources which we assert are most 
likely to be classical Be stars, "Type-2:ES" sources (green triangles) correspond to our polarimetrically designated Type-2 sources which 
exhibit electron scattering intrinsic polarimetric signatures and are thus likely classical Be stars, "Type- 2-4, other" sources (dark blue squares) 
correspond to all other candidate Be stars in clusters investigated in this paper regard less of whether they were detected via our polarimetric 
survey or not, "ELHC Be" sources (yellow triangles) correspond to LMC stars which Ide Wit et al.l 1120 051) claim are likely to be classical Be 
stars, "Galactic Be" sources (black triangles) are 101 known Galactic Be st ars tabulated by Dougherty ct al. ( 199l]) , "MC B[e]" (light blue 
cross) correspond to 2-MASS and ground-based (Gummcrsbach ct al. 1995) photometry of known Magellanic Cloud B[e] supergiants, and 
"MC HAeBe" (pink triangle) corresponds to the star ELHC-7. which Ide Wit et all J2005D assert is a Herbig Ae/Be star which resides in the 
LMC. As discussed in Section 15.3.31 most of the candidate Be stars which we were able to correlate with the 2MASS catalog exhibit colors 
consistent with those expected from classical Be stars, and inconsistent with those expected from composite gas plus dust disk systems (i.e. 
Herbig Ae/Be or B[e] stars). 
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Fig. 16. — The same data as presented in Figure [T5l without error bars. 
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Fig. 17. — The subset of data presented in Figure ITrJl corresponding to candidate Be stars residing in "very young" clusters of age 5-8 
Myr, is presented. As discussed in Section 15.3.31 the data presented in this Figure and Figure [TBI exhibit similar properties. Regardless of 
polarimetric classification type, all candidate Be stars in clusters of age 5-8 Myr appear consistent with the expected colors of classical Be 
stars and inconsistent with the observed colors of composite gas plus dust systems (i.e. Herbig Ae/Be, B[e] stars). We interpret these data 
as additional evidence that many of the photometrically identified candidate Be stars in clusters 5-8 Myr old by WB06 are classical Be stars 
and not pre-main-sequence Herbig Be stars. 



